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1 Introduction 

The standard model (SM) of particle physics has been extremely successful in describing phe- 
nomena at the highest energies attained thus far. Nevertheless, it is widely believed to be only 
an effective description of a more complete theory, which supersedes it at higher energy scales. 
Of particular theoretical interest is supersymmetry (SUSY) [1-6], which solves the "hierarchy 
problem" (ZUSI of the SM at the expense of introducing a large number of supersymmetric par- 
ticles with the same quantum numbers as the SM particles, but differing by half a unit of spin. 
If R-parity conservation [9J is assumed, supersymmetric particles are produced in pairs and de- 
cay to the lightest supersymmetric particle (LSP). If the LSP is neutral and weakly-interacting, 
it goes undetected giving rise to a signature with missing energy. 

Experiments at the energy frontier, i.e. at the Fermilab Tevatron collider ||T0l - [T3| and previously 
at the CERN SppS HI [^, HERA HH HZI and LEP ^ colliders, have performed extensive 
searches for signs of SUSY. In the absence of a positive signal, lower limits on the masses of 
SUSY particles have been set. With its higher centre-of-mass energy of 7 TeV, the Large Hadron 
Collider (LHC) at CERN could produce SUSY particles (sparticles) with masses larger than 
the current limits. The dominant production channels of heavy coloured sparticles at the LHC 
are squark-squark, squark-gluino and gluino-gluino pair production. In the context of SUSY 
with R-parity conservation, heavy squarks and gluinos decay into quarks, gluons and other 
SM particles, as well as a neutralino (i.e. the LSP), which escapes undetected, leading to final 
states with several hadronic jets and large missing transverse energy. While squark-squark 
production usually leads to two jets, gluino production typically results in higher jet multiplic- 
ities. This Letter describes a search for the production and decay of SUSY particles by the CMS 
experiment, in events with two or more energetic jets and significant imbalance of transverse 
energy. 

The search is not optimized in the context of any particular model of SUSY. To interpret the 
results, a simplified and practical model of SUSY-breaking, the constrained minimal supersym- 
metric extension of the standard model (CMSSM) ['T9'.'20], is used. The CMSSM is described by 
five parameters: the universal scalar and gaugino mass parameters (nio and mi/2, respectively), 
the universal trilinear soft SUSY breaking parameter Aq, and two low-energy parameters, the 
ratio of the two vacuum expectation values of the two Higgs doublets, tan fi, and the sign of the 
Higgs mixing parameter, sign(/^). Throughout the Letter, two CMSSM parameter sets, referred 
to as LMO and LMl [21 J, are used to illustrate possible CMSSM yields. The parameter values 
defining LMO are mo = 200GeV,?«i/2 = 160GeV, Aq = -400GeV,tan^ = 10, and sign (/^O > 0. 
Those for LMl are mo = 60 GeV, mi/2 = 250 GeV, Aq = 0, tan j6 = 10, and sign(ju) > 0. 

2 The CMS Detector 

The central feature of the CMS apparatus is a superconducting solenoid, 13 m in length and 
6 m in diameter, which provides an axial magnetic field of 3.8 T. The bore of the solenoid is in- 
strumented with various particle detection systems. The steel return yoke outside the solenoid 
is in turn instrumented with gas detectors used to identify muons. Charged particle trajecto- 
ries are measured by the silicon pixel and strip tracker, with full azimuthal coverage within 
\rj\ < 2.5, where the pseudorapidity f/ is defined as t] = — lntan(0/2), with 6 being the polar 
angle of the trajectory of the particle with respect to the counterclockwise beam direction. A 
lead-tungstate crystal electromagnetic calorimeter (ECAL) and a brass /scintillator hadron cal- 
orimeter (HCAL) surround the tracking volume and cover the region \ri\ < 3. In the region 
< 1.74, the HCAL cells have widths of 0.087 in pseudorapidity and 0.087 in azimuth {(p). In 
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the {rj, (p) plane, and for \ri\ < 1.48, the HCAL cells map on to 5 x 5 ECAL crystal arrays to form 
calorimeter towers projecting radially outwards from close to the nominal interaction point. At 
larger values of the size of the towers increases and the matching ECAL arrays contain 
fewer crystals. Within each tower, the energy deposits in ECAL and HCAL cells are summed 
to define the calorimeter tower energies, subsequently used to provide the energies and di- 
rections of hadronic jets. The detector is nearly hermetic, which allows for energy-balance 
measurements in the plane transverse to the beam axis. A more detailed description of the 
CMS detector can be found elsewhere p2|. 

3 Event Selection 

3.1 Hadronic final state selection 

The data sample used in this analysis is recorded with a trigger based on the scalar sum of the 
transverse energy Ej of jets, defined in general as Hj = J^J^l Ej^', where Njet is the number of 
jets. Events are selected if they satisfy Hj^'^^^^^^ > 150 GeV. 

At the trigger level, the calorimeter response is not corrected to achieve a uniform and absolute 
scale of transverse jet energy; nevertheless the trigger requirement is fully efficient for events 
with an offline-reconstructed Hj in excess of 250 GeV, thus providing a high signal efficiency 
for the region of the CMSSM parameter space relevant for the present search, where squarks 
and gluinos have masses of several hundred GeV. Additionally, events are required to have at 
least one good reconstructed pp interaction vertex Il23ll . 

Jets are reconstructed offline from the energy deposits in the calorimeter towers, clustered by 
the anti-kx algorithm Il24l with a size parameter of 0.5. In this process, the contribution from 
each calorimeter tower is assigned a momentum, the magnitude and direction of which are 
given by the energy measured in the tower and the coordinates of the tower. The raw jet energy 
is obtained from the sum of the tower energies, and the raw jet momentum by the vectorial sum 
of the tower momenta, resulting in a nonzero jet mass. The raw jet energies are corrected to 
establish a relative uniform response of the calorimeter in rj and a calibrated absolute response 
in transverse momentum pj. The uncertainty on the energy scale of these corrected jets varies 
between 3% and 5%, depending on the jet pj and \r]\ |25| . The jets considered in this analysis 
are required to have Ej > 50 GeV, | // 1 < 3 and to pass jet identification criteria 126} designed 
to reject spurious signals in the calorimeters. The pseudorapidity of the jet with the highest 
Ej (leading jet) is required to be within \rj\ < 2.5 and the transverse energy of each of the two 
leading jets must exceed 100 GeV. 

Events with jets passing the Ej threshold but not satisfying the jet identification criteria or 
the rj acceptance requirement are vetoed, as this deposited energy is not accounted for in the 
event kinematics. Similarly, events in which an isolated lepton (electron p7) or muon 1128 1 ) 
with Pj > 10 GeV is identified are rejected to suppress events with genuine missing energy 
from neutrinos. Furthermore, to select a pure multi-jet topology, events are vetoed in which 
an isolated photon p9| with pj > 25 GeV is found. These vetoes reject 5% of the previously 
selected events in data and simulation. 

At this preselection stage, the background from multijet production, as predicted by quantum 
chromodynamics (QCD), is still several orders of magnitude larger than the typical signal ex- 
pected from SUSY. The Hj distribution for the selected events is shown in Fig.jljand compared 
to simulation-based background estimates. The QCD multijet background is estimated using 
the PYTHIA 6.4 ||30| Monte Carlo generator with tune Z2 [31 1 . Electro weak backgrounds from 



3.2 Final event selection for SUSY search 



3 




500 1000 1500 2000 2500 

K,(GeV) 



Figure 1: Hj distribution after preselectior\, for data as well as for all star\dard model back- 
grounds and two SUSY signal samples with parameter sets LMO and LMl, normalized to an 
integrated luminosity of 35 pb^^. The hatched area corresponds to the uncertainty in the SM 



estimate as defined in Section 3.1 The SM distributions are only displayed for illustration 
purposes, as they are the result of Monte Carlo simulation, while the actual estimate of the 
background from SM processes in this search is based on data, as described in detail in Sec- 
tion 141 



W + jets, Z — )■ 1/1/ + jets and tt + jets events, which will be referred to collectively as the elec- 
troweak (EWK) backgrounds in what follows, are simulated using MadGraph [32|. The SM 
distribution, i.e. the sum of the QCD multijet and EWK distributions, is indicated in Fig. [T] 
as a hatched band representing the combined statistical and systematic uncertainties from the 
jet energy scale and resolution. The expected Hj distributions for two low-mass SUSY signal 
points, LMO and LMl, are overlaid. With the exception of tt, the SM processes fall off expo- 
nentially over the entire Hj range, whereas a broad peak at values of a few hundred of GeV 
is expected for the signal models. The selection is tightened by requiring the Hj of all jets to 
exceed 350 GeV, thus ensuring large hadronic activity in the event. This requirement substan- 
tially reduces the contributions from SM processes while maintaining a high efficiency for the 
SUSY topologies considered. 

3.2 Final event selection for SUSY search 

Jet mismeasurements, caused by possible detection inefficiencies or by nonuniformities in the 
calibration of the calorimeters, are the dominant source of large missing transverse energy £ j 
in events from QCD multijet production. To control this background and to separate it from a 
genuine missing energy signal, a variable that is robust against energy mismeasurements, ocj, 
is used. For events with two jets, uj, first introduced in Refs. I|2TI|33| and inspired by Ref. ||34| , 
is defined as 



where Ej'^ is the transverse energy of the less energetic of the two jets in the event and Mj is 
the transverse mass of the di-jet system, defined as 



4 



3 Event Selection 



For a perfectly measured di-jet event, with Ej^^ = £t'^ and jets back to back in (p, and in the 
limit where the jet momenta are large compared to their masses, the value of aj is 0.5. In the 
case of an imbalance in the measured transverse energies of back to back jets, ccj takes on values 
smaller than 0.5, while for jets that are not back to back, ocj can be greater than 0.5. 

For larger jet multiplicities, the n-jet system is reduced to a di-jet system by combining the jets 
in the event into two pseudo-jets. The Ej of each of the two pseudo-jets is calculated as the 
scalar sum of the contributing jet Ex's. The combination chosen is the one that minimizes the 
Ej difference between the two pseudo-jets. This simple clustering criterion has been found to 
result in the best separation between QCD multijet events and events with genuine Bj- 

Values of aj above 0.5 can occur for QCD multijet events, either with multiple jets failing the 
Ej > 50 GeV requirement, or with missing transverse energy arising from jet energy resolution 
or severe jet energy under-measurements due to detector inefficiencies. On the other hand, 
events with genuine £ j often have much larger values of ocj, resulting in a good separation of 
signal events from the QCD multijet background. 

The Kj distributions are shown separately for di-jet and > 3-jet events in Fig. |2j As antici- 
pated, these distributions peak at ocj = 0.5 for QCD multijet events and then fall sharply in the 
range 0.5 to 0.55, reaching a level 4-5 orders of magnitude lower than the peak value. Multijet 
events from QCD background are therefore efficiently rejected by requiring ocj to exceed 0.55. 
Given the selection requirement Hj > 350 GeV, this threshold on kj is equivalent to demand- 
ing Ht/Ht > 0.4, i.e.to Hi > 140 GeV. 

To reject events with false missing energy arising from significant jet mismeasurements in 
masked regions of the ECAL, which amount to about 1% of the ECAL channel count, the 
following procedure is employed. The jet-based estimate of the missing transverse energy, 
Ht = Ifljl = I — Ejets ?^jetl' which is obtained by summing the transverse momenta of all the 
jets in the event, is now recomputed while ignoring one of the reconstructed jets. The differ- 
ence in azimuth between the recomputed Hi and the ignored jet is then calculated. The fij is 
recomputed for each configuration that results from ignoring, in turn, each of the jets in the 
event, while leaving all other jets intact, and the minimum of all the azimuthal differences, 
A^*, is found. The jet whose subtraction from the calculation yields this minimum value, 
is identified as the jet that is most likely to have given rise to the in the event. Events with 
A^* < 0.5 are rejected if the distance in the {rj, (p) plane between the selected jet and the closest 
masked ECAL region, ARecal/ is smaller than 0.3. 

Artificially large values of Hj can also result in events with multiple jets below the selection 
requirement of Ej > 50 GeV, since these jets are not included in the computation of Hi- To 
protect against these events, flj, i.e. the jet-based estimate of the missing energy. Hi, is com- 
pared to the calorimeter tower-based estimate, £ which includes the energy from all jets, 
irrespective of threshold p5| . Events with Rmiss = Bj/lj^^° > 1.25 are rejected. 

Table [l] lists the number of events passing each step of the event selection for data and simu- 
lation. The expectations from simulation are listed only for comparison; the actual expected 
yields from standard model processes are determined from control data samples, as described 
in the following section. 
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Figure 2: Distribution of aj for di-jet events (left) and > 3-jet events (right), requiring Hj > 
350 GeV. Events with (Xj > 1.5 are included in the rightmost bin. In both figures the hatched 



area corresponds to the uncertainty in the SM estimate as defined in Section 3.1 



After the selection requirements on (Xj, ARecal arid Kmiss/ the QCD multijet background pre- 
dicted by PYTHIA 6.4 is less than one event for an integrated luminosity of 35 pb^^. This esti- 
mate is also obtained with PYTHIA 8.1 [36 1 (tune 1) and with the MadGraph generator. After 
all selection requirements, the only significant remaining background stems from electroweak 
processes with genuine Bjiri the final state. In the di-jet case, the largest backgrounds with real 
missing energy are the associated production of W or Z bosons with jets, followed by the weak 
decays Z — > t/t7 and W — >• tv, or by leptonic W/ Z decays in which one or more leptons are not 
reconstructed. At higher jet multiplicities, tt production followed by semileptonic weak decays 
of the t and t quarks becomes important. In this case, the three backgrounds, Z — > i/t7 + jets, W 
+ jets and tt, are of roughly equal size. The largest fraction of the W + jets and tt backgrounds 
stem from W — )■ tv decays where in two thirds of the cases the t decays hadronically and is 
identified as a jet. The two remaining backgrounds from electrons or muons produced in W 
decays that fail either the isolation or acceptance requirements (px > 10 GeV and t] coverage) 
are of similar size. 



4 Background Estimate from Data 

The SM background in the signal region is estimated directly from data using two independent 
methods. The first method makes use of control regions at lower Hj to estimate the total back- 
ground from all SM processes (Section [4.1| , while the second method estimates the contribution 
from electroweak processes using W — ?► + jets (Section [42] | and 7 + jets (Section |43] l events 
in the data. 

4.1 Inclusive background estimate 

The total background can be estimated from two control regions at low Hj: the HT250 region, 
which contains events with Hj between 250 and 300 GeV, and the HT300 region, which contains 
events with Hj between 300 and 350 GeV. Given the current experimental limits on the squark 
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Table 1: The number of events observed and expected from Monte Carlo simulation after the se- 
lection requirements, for data and background samples (QCD multijet simulated with PYTHIA 
6.4(Z2), Z — > vv, W +jets, tt). The quoted errors represent the statistical uncertainties on the 
yields and all numbers are normalized to an integrated luminosity of 35 pb^^. 



Selection 


Data 


SM 


QCD multijet 


Z — > vv 


W+jets 


tt 


Ht > 250 GeV 


4.68M 


5.81M 


5.81M 


290 


2.0k 


2.5k 


Et'" > 100 GeV 


2.89M 


3.40M 


3.40M 


160 


610 


830 


Ht > 350 GeV 


908k 


l.llM 


l.llM 


80 


280 


650 


OCT > 0.55 


37 


30.5±4.7 


19.5±4.6 


4.2±0.6 


3.9±0.7 


2.8±0.1 


AKecal > 0.3 V A^* > 0.5 


32 


24.5±4.2 


14.3±4.1 


4.2±0.6 


3.6±0.6 


2.4±0.1 


^miss ^ 1.25 


13 


9.3±0.9 


0.03±0.02 


4.1±0.6 


3.3±0.6 


1.8±0.1 



and gluino masses, these two regions are expected to be dominated by SM processes. The 
search region for the signal, which is referred to as the HT350 region in what follows, is defined 
as events with Hj > 350 GeV. 

The method is based on the variable Raj, defined as the ratio of the number of events passing 
and failing a requirement on ccj, given all other selection requirements. To minimize the effi- 
ciency bias arising from the phase space reduction in the lower Hj regions, the pj thresholds 
for the two bins are adjusted to keep the ratio of pj/ Hj constant in each region. In the HT300 
region, the resulting thresholds are 86 GeV for the two leading jets and 43 GeV for additional 
jets. In the HT250 region the respective thresholds are 71 and 36 GeV. In the absence of a SUSY 
signal, the ratio R^t can then be extrapolated from the measured values in both control regions 
to predict the value in the signal region, HT350. 

Figure |3] (left) shows the evolution of Raj as a function of Hj for two thresholds on (Xj, namely 
a.j > 0.51 and ocj > 0.55. For aj > 0.51, the numerator of Raj is dominated by the QCD 
multijet background, for which the missing transverse energy mostly originates from energy 
mismeasurements. As the relative resolution of calorimetric energy measurements improves 
with energy, and therefore with Hj, the relative importance of this background is expected to 
decrease with increasing Hj. This effect is clearly visible in Fig. |3] (left), which shows the falling 
behaviour for seven equidistant bins in Hj. In contrast, for ax > 0.55, the numerator of Raj is 
dominated by the electroweak background, with genuine ^ t, the relative importance of which 
is expected to be constant with increasing Hj. 

The latter behaviour is confirmed in an independent sample of events with a W decaying to 



}iv, accompanied by jets. (The selection of these events is given in Section 4.2 ) The ratio of 
the number of selected W — > /^(v + jets events to the number of events failing the W selection 
(hence dominated by the QCD multijet background) is shown in Fig. |3] (left) for the same Hj 
bins, and confirms the independence of Raj on Hj when the numerator is dominated by events 
with genuine ^j. 

The ratio Raj in the HT350 region can be estimated from the Raj values measured in regions 
HT250 and HT300, for ccj > 0.55, using the double ratio Rr: 

^ R,,(HT300) ^ R,,(HT350) 
^ Ra^(HT250) Ra^(HT300)' ^ ' 
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Figure 3: Evolution of the ratio Raj as a function of Hj for events with Njet > 2; (left) for data 
and SM backgrounds, and two different values of ocj, as well as for an independent W — > f/i/ + 
jets control sample (Section |4.2| ; (right) for the SM backgrounds added to the SUSY signal ex- 
pected from each of the two benchmark points, LMO and LMl. Markers are offset horizontally 
for improved visibility. 

The total number of events with ocj > 0.55 expected from SM processes in the signal region 
is the product of the extrapolated Rfl;T(HT350) and the number of events with ocj < 0.55 in 
the Hj > 350 GeV region. The total number of background events in HT350 thus estimated is 
9-4^4 Q (stat) ±1.0 (syst). The dominant systematic uncertainty for this method is estimated by 
varying the relative magnitude of the three EWK processes, while maintaining the Hj depen- 
dence of each process to the one predicted by the Monte Carlo. The background estimate is 
insensitive to this variation: since Rr ~ 1, the change in the estimate is always much smaller 
than the statistical uncertainty. Even under extreme variations of the individual EWK processes 
by up to five times higher values than those predicted by simulation, the systematic uncertainty 
is at least a factor two smaller than the statistical uncertainty. The same comments hold for the 
background estimate variants described below. 



For comparison, Fig.|3](right) shows the expectation of Raj if a SUSY signal from each of the two 
benchmark points LMO and LMl were present in addition to the SM backgrounds. The signal 
is predominantly visible in the HT350 region, especially for LMl, which has larger squark and 
gluino masses. There is, nevertheless, a sizeable signal in the HT300 control region as well. 
This contamination of the control region by a potential signal is taken into account in the limit 
calculation (Section |5.1| an d both benchmark points are ruled out with a confidence level of 



99% or higher (Section 5.2 1. 



A variant of this background estimation method relies on the Raj measurement in the HT300 
bin only and uses a small correction from MC simulation to predict R^j in the signal region. 
This variant results in an estimate of 12.0±g 3 (stat) ±0.4 (syst) events. Another variant, also 
based on the independence of Raj on Hj when the data sample is dominated by EWK pro- 
cesses, i.e. for OLj > 0.55, uses the weighted average of the Raj values measured in the two 
control regions. This value is then also used in the signal region to obtain a background esti- 
mate of 12.5±1.9 (stat) ±0.7 (syst). Within uncertainties, the three estimates are in agreement. 
Furthermore, in the simulation all methods are shown to provide an unbiased estimator of 
the number of total background events. For the remainder of this Letter the result of the first 
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method, which relies entirely on measured data and makes the most conservative assumption 
on the evolution of the double ratio with Hj, is used to estimate the total background. 

4.2 W + jets and tt background 

A second background estimation method uses an independent selection of W — ?► //v + jets 
events in the data in order to assess the contribution from SM processes with genuine Kt- The 
W ^ jiv + jets are selected as described in [37], with an energetic and isolated muon in the final 
state, and by requiring the transverse mass of the W to be larger than 30 GeV (to ensure a very 
pure sample originating from W + jets and tf). The muons are required to be separated from 
the jets in the event by a distance larger than 0.5 in the (?/, (p) plane. Since aj > 0.55 implies 
Hj/ Hj > 0.4, only events with Ht > 140 GeV are considered in the signal region (HT350). In 
the lower Hj regions, this requirement is scaled accordingly to Ht > 120 (100) GeV for HT300 
(HT250). 

In the HT350 region this selection yields 25 events, in agreement with the 29.4 ±1.4 events 
predicted by the simulation. In the HT250 (HT300) region, 134 (52) W candidates are recon- 
structed, in agreement with the prediction of 135.5 ± 3.2 (56.7 ± 2.2) events. The fraction of 
W ^ j-iv + jets events with ocj > 0.55 in the data is also in good agreement with the simulation: 
seven data events are found in the signal region, compared with 5.9 ± 0.6 events predicted, 
whereas 32 (12) events in the data pass the ocj > 0.55 requirement in the HT250 (HT300) region, 
compared to 29.2 ± 1.4 (11.1 ± 1.1) events expected. 

The number of W + jets and tt events satisfying the hadronic final state selection of Section |3| 

^dl'ta^^'^' estimated from the number of events in the muon sample, N^'J', and the ex- 

pected relative ratio of these two types of events. The value of this ratio is taken from Monte 
Carlo simulation, which yields N^;,^^'^ = N^^^""^ / N^^^' x N^;^^ ^ 0.86 x N^;f . The total 
background from W + jets and tt processes is thus estimated to be 6.ll^| (stat) ± 1.8 (syst). 

Given the reliance on simulation for the factor N^q^^'^ / N^q" , conservative uncertainties on all 
the parameters entering this ratio have been assigned. The systematic uncertainty is estimated 
to be 30% and is dominated by the uncertainty on the efficiency for vetoing leptons. 

4.3 Z VV+ jets background 

The remaining irreducible background stems from Z — ?► i/i7 + jets events. An estimate of this 
background can be obtained from 7 + jets events, which have a larger production cross sec- 
tion but kinematic properties similar to those of Z — ?► vv + jets events when the photon is 
ignored [38j. These 7 + jets events provide a measurement of the acceptance of the uj > 0.55 
requirement directly from data. The 7 + jets sample is selected by requiring photons, i.e. lo- 
calized electromagnetic depositions satisfying very tight isolation criteria, with pj > 100 GeV, 
< 1.45, and with a distance in the {tj,(p) plane to any jet larger than 1.0. Subsequently, the 
photon is ignored and the same hadronic final state selection as described in Section |3] is ap- 



plied. As in Section 4.2 Ht is required to exceed 140 GeV. This selection yields seven events in 
the data compared with 6.5 ± 0.4 expected from simulation. The relative acceptances, together 
with the appropriate ratio of cross sections for 7 + jets and Z — >• vv + jets, taken from simula- 
tion, are then used to estimate the number of Z — > vv + jets events in the signal region, found 
to be N{Z — >■ vv + jets) = 4.41^^ (stat) ± 1.8(syst). The main systematic uncertainties arise 
from the ratio of cross sections between 7 + jets and Z — > viz + jets in the simulation (30%), the 
efficiency for photon identification (20%), and the purity of the photon selection (20%), which 
add up to ^ 40%. 



4.4 Background estimate summary 
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To check the validity of this uncertainty estimate, the number of 7 + jets events can also be used 
to predict the number of W + jets events. Ten W + jets events are observed with 250 < Hj < 
350 GeV, Njets = 2 and ocj > 0.55, in agreement with the prediction of 8.5 ± 1.5(stat) ± 2.6(syst) 
from the 7 + jets process. This agreement gives cor\fidence that the magnitude of the assigned 
systematic uncertainties is adequate. 

As a further cross-check, the W — > + jets sample discussed above is used to estimate the 
background from Z — > 1/1/ + jets events. With the observed number of reconstructed W — > f/v 
+ jets events, the ratio of cross sections and branching ratios for W and Z bosons, and the ratio 
of reconstruction efficiencies estimated from simulation, 4.9^^ g (stat) ±1.5 (syst) Z — > i/i7 + jets 
events are predicted, in agreement with the value obtained from the 7 + jets sample. 

4.4 Background estimate summary 

The SM backgrounds to this analysis have been evaluated with independent data control sam- 
ples. First, from the lower Hj regions in the data, a prediction for the total SM background 
of 9.4^^ Q (stat) ±1.0 (syst) events in the signal region is obtained. Then, with a W — > jiv + 
jet control sample, a contribution of 6.1^j g (stat) ±1.8 (syst) events from the combination of 
W + jets and tt processes is estimated. Finally, with the 7 + jets sample, the background from 
Z —> vv + jets events is estimated to be 4.4^^ ^ (stat) ±1.8 (syst). Therefore, the estimate of the 
SM background arising from EWK processes with genuine £t is 10.5^2 5 events, which is in 
agreement with the inclusive estimate obtained from the lower Hj control regions. All these 
background estimates are used in the limit calculation. 

The potential effect of multiple interactions per biinch crossing (pileup) in the event selec- 
tion and in the background estimates is evaluated by comparing the fraction Ra^, for several 
thresholds on ccj, between events with only one primary vertex and events with more than one 
primary vertex. No discernible difference has been found. 

The remaining events are found to exhibit the topological and kinematic properties expected 
from the SM backgrounds. Two distributions, which are expected to show good separation 
between the SM background and the SUSY signal, are shown in Fig. |4] The A^* variable (Sec- 
tion is useful in identifying mismeasured jets, since jet mismeasurements in QCD multijet 
events result in small values of A^*, whereas events with genuine ^j, e.g. from EWK processes, 
populate A^* evenly. In Fig. |4] (left) the A^* distribution for the 13 data events which pass all 
selection requirements is displayed. The data are consistent with EWK processes and there is 
no indication of an enhanced contribution from QCD multijet processes which would manifest 
itself at small values of A^* . 

The "effective mass", Mgff = Hj ± ^^j, which characterizes the overall energy scale of the event, 
is shown in Fig. |4] (right) after all selection requirements. The data are compared with the SM 
background expectation along with two SUSY benchmark points. The shape and magnitude 
of the Meff distribution observed in the data are consistent with the expectation from the SM 
backgrounds. The yields expected from the LMO and LMl benchmark SUSY models are in 
excess of the data over most of the Mgff range. 

Both these variables exhibit differences between SUSY signal events and events from SM back- 
grounds and could, therefore, be used to improve the limits extracted in the following section. 
We have chosen not to do so because the current search has been optimized for the demonstra- 
tion of a potential new signal, rather than for the extraction of the most stringent limits in the 
SUSY parameter space. 

In summary, 13 events are observed in the data, a yield consistent, within the uncertainties. 
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5 Interpretation of the Result 



with the expectation from the SM processes. In addition, the kinematic properties of these 
events are consistent with the EWK backgrounds, with a negligible contribution from QCD 
multijet processes. 




Figure 4: Left: the Ai^* distribution after all selection requirements. Right: the effective mass 
Meff distribution after all selection requirements for SM processes and two low mass SUSY 
benchmark points. In both figures the hatched area corresponds to the uncertainty in the SM 



estimate as defined in Section 3.1 



5 Interpretation of the Result 
5.1 Method and limit on signal yield 

The background estimation methods described in the previous section are combined to provide 
an estimate of the total number of background events. This estimate is found to be compatible 
with the number of events selected. An upper limit on the number of non-SM events con- 
sistent with the measurements is derived using the Feldman-Cousins method [39], which is 
generalized to take into account nuisance parameters by using the Profile Likelihood ratio flO]. 
The input to the profiling method is the total likelihood function L^Q^^l for the measurements 
in the control and signal regions. To construct this likelihood function, the numbers of events 
observed in the signal region and in each of the control samples are treated as independent 
event-counting experiments. The total likelihood function can be written as 



^total — ^signal ' ^inclusive ' ^W/tt ' ^Z^vvr (2) 

where the different factors correspond to the likelihood of the measurement in the signal re- 
gion, Lsignai, of the inclusive background extraction method using all three Hj bins defined 



in Section 4.1 Lin^iusive/ of the exclusive measurement of tt and W + jet event background de- 



scribed in Section 4.2 L^j /ti, arid of the Z ^ vv background component described in Section |43} 

The likelihood functions are taken as Poisson-distributed probabilities to measure the number 
of events observed while expecting to see the estimated number of background events plus 
a certain fraction of CMSSM SUSY signal events. The expected backgrounds in the auxiliary 



5.2 Interpretation within tine CMSSM 
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measurements and in the signal-like region are related as described in Section |4j The ratio of 
the expected signal events in the control region and in the signal region is model-dependent 
and varies from point to point in the CMSSM parameter space. 

The expected number of background events and the systematic uncertainties on the back- 
ground prediction and on the signal selection efficiency are treated as nuisance parameters. 
The probability density functions which describe the systematic uncertainties are assumed to 
be Gaussian with variance given by the systematic uncertainties derived in the previous sec- 
tion. 

The systematic uncertainties on the signal event yield can be split into two parts: theoretical un- 
certainties on the predicted cross section of the different production processes (squark-squark, 
squark-gluino, gluino-gluino) and experimental uncertainties on the integrated luminosity [41] 
and on the selection efficiency. 

Whereas the theoretical uncertainties are strongly model dependent, the experimental uncer- 
tainties are found to be essentially independent of the signal model. The experimental sys- 
tematic uncertainties on the estimated event yield are the uncertainty on the luminosity mea- 
surement (11%), the effect of rejecting events with jets pointing to masked ECAL regions (3%), 
the modelling of the lepton and photon vetoes in the simulation (2.5%), and the effect of the 
uncertainty in the jet energy scale and resolution on the selection efficiency (2.5%). These un- 
certainties are included in the limit calculation. The effect of multiple interactions per bunch 
crossing on the signal is evaluated by comparing the efficiency for signal events passing all se- 
lection requirements with and without the inclusion of multiple interactions in the simulation. 
The effect on the efficiency is negligible. 

If a potential signal contamination in the background control samples is ignored, an upper 
limit on the number of signal events compatible with the observations at 95% confidence level 
(CL) can be obtained. For an integrated luminosity of 35 pb^^ this number is 13.4 events. The 
p value for the hypothesis of standard model background only, calculated from the ratio of 
likelihoods, is 0.3. 

5.2 Interpretation within the CMSSM 

To interpret the consistency of the observed number of events with the background expectation 
in the context of a model, and also to facilitate the comparison with previous experimental 
results, an exclusion limit in the CMSSM is set. This limit is obtained by testing, for each 
point in the parameter space, whether the number of signal events predicted after all selection 
requirements is compatible with observations at 95% CL. 

Signal contamination in the data control samples used to estimate the background is also taken 
into account by explicitly including the number of signal expected in the control regions. As 
the search is designed for robustness and background control, the same selection is applied at 
each point in the parameter space, and no dedicated optimization for the CMSSM parameter 
space is performed. 

An example of the analysis efficiency and corresponding event yields after all selection require- 
ments, broken down by the most relevant production processes (squark-squark (qq), squark- 
gluino (qg), and gluino-gluino (gg)), is presented in Table |2]for the benchmark point LMl. Two 
different experimental efficiencies etotal arid ^signature are given. The first number, ^total/ is nor- 
malized to the total number of signal events in LMl, while esignature is defined with respect 
to the total number of all-hadronic events in LMl where, as in the analysis, leptons and pho- 
tons are vetoed. For the different production mechanisms, ^totai varies from 12% to 16%. The 
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6 Summary 



Table 2: Breakdown of expected event yields and selection efficiencies for the most important 
production channels of the LMl benchmark point after all selection requirements. No distinc- 
tion has been made between q and q. The quoted errors represent the statistical uncertainties 
on the yields and efficiencies. The efficiencies etotai srid ^signature ^re defined in the text in Sec- 
tion |52] 



Production mechanism 


Yields for 35 pb^^ 


£total(%) 


^signature ( ^) 


qq 


9.7±0.1 


16.0±0.1 


22.2±0.4 


qg 


8.8±0.1 


14.4±0.1 


23.0±0.5 


gg 


0.71±0.02 


12.0±0.4 


22.5±2.0 



signature-based efficiency is almost constant, varying between 22% and 23%, which indicates 
that the analysis has a uniform sensitivity to the different production channels in LMl. With 
the current data, the LMl and LMO benchmark points are excluded at 99.2% CL and 99.99% 
CL, respectively. 

The 95% CL limit in the (mo,mi/2) plane, for tan j6 = 3, Aq = and sign(f/) > 0, is shown in 
Fig. |5j The SUSY particle spectrum is calculated using SoftSUSY [42J, and the signal events are 
generated at leading order (LO) with pythia 6.4. Next-to-leading order (NLO) cross sections, 
obtained with the program Prospino Il43l , are used to calculate the observed and expected 
exclusion contours. Systematic uncertainties on the NLO predictions due to the choice of the 
renormalization and factorization scales have been taken into account. The uncertainties on 
the used parton distribution functions (PDF) for CTEQ6.6 [44] are estimated from the envelope 
provided by the CTEQ6.6 error function. For reference, the observed limit using LO cross 
sections is also shown. 

The expected limit covers a larger part of the [nio, mi/2) plane than the measured limit, as the 
number of events observed in the signal region is slightly larger than the number of background 
events predicted from the control regions. The excluded regions for the CDF search for jets + 
missing energy final states [TOl were obtained for tan j5 = 5, while those from DO [12] were 
obtained for tan /3 = 3, each with approximately 2 fb^^ of data. The LEP-excluded regions 
are based on searches for sleptons and charginos [18J. A comparison of the exclusion limit for 
tan /3 = 3 to that for tan j6 = 10 for fixed values of Aq = and sign(/^) > indicates that the 
exclusion reach is only weakly dependent on the value of tan j6; the limit shifts by less than 
20GeV in mo and by less than lOGeV in mi/2. The DO exclusion limit, valid for tan/5 = 3 
and obtained from a search for associated production of charginos ^rid neutralinos Xi 
trilepton final states 1, 13 J , is also included in Fig. |5j In contrast to the other limits presented 
in Fig. |5| the result of the trilepton search is strongly dependent on the choice of tan and it 
reaches its highest sensitivity in the CMSSM for tan jS values below ten. 

6 Summary 

The first search for supersymmetry in events collected by the CMS experiment from proton- 
proton collisions at a centre-of-mass energy of 7 TeV has been presented. The final states with 
two or more hadronic jets and significant missing transverse energy, as expected from high- 
mass squark and gluino production and decays, have been analysed in data corresponding to 
an integrated luminosity of 35 pb^^. A search for a SUSY signal has been performed at high val- 
ues of the scalar sum of the transverse energy of jets, Hj. The primary background, from QCD 
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100 200 300 400 500 

mo (GeV) 

Figure 5: Measured (red line) and expected (dashed blue line) 95% CL exclusion contour at 
NLO in the CMSSM (mo,mi/2) plane for tan/3 = 3, Aq = and sign(^) > 0. The measured 
LO exclusion contour is shown as well (dot-dashed green line). The area below the curves is 
excluded by this measurement. Exclusion limits obtained from previous experiments are pre- 
sented as filled areas in the plot. Grey lines correspond to constant squark and gluino masses. 
The plot also shows the two benchmark points LMO and LMl for comparison. 

multijet events, has been reduced by several orders of magnitude down to a negligible level 
using a robust set of requirements designed specifically for the exploratory early data-taking 
phase of the experiment. The sum of standard model backgrounds has been estimated from an 
extrapolation of the data observed at lower Hj values. The only remaining backgrounds have 
been found to stem from electroweak processes, namely W + jet, Z + jet, and tt production, 
where the weak decays of the vector bosons involve high-momentum neutrinos. An indepen- 
dent estimate of the electroweak backgrounds, from W — ?► f/t/ + jets decays as well as 7 + jets 
events in the data together with input from simulation, has been found to be well compatible 
with the estimate from control samples in data. Here, conservatively large systematic uncer- 
tainties have been assigned to the background estimates. The measurements are in agreement 
with the expected contributions from standard model processes. Limits on the CMSSM param- 
eters have been derived, and have been shown to improve significantly those set by previous 
experiments. 
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